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ABSTRACT 

We stack Spitzer 24 /im images for ~ 7000 galaxies with O.I < z < I in the Chandra Deep Field 
South to probe the thermal dust emission in low-luminosity galaxies over this redshift range. Through 
stacking, we can detect mean 24 /rm fluxes that are more than an order of magnitude below the 
individual detection limit. We find that the correlations for low and moderate luminosity galaxies 
between the average Lir/Luv and rest-frame i?-band luminosity, and between the star formation 
rate (SFR) and Lir/LuV) are similar to those in the local Universe. This verihes that oft-used 
assumption in deep UV/optical surveys that the dust obscuration-SFR relation for galaxies with SFR 
< 2OM0yr“^ varies little with epoch. We have used this relation to derive the cosmic IR luminosity 
density from 2; = 1 to z = 0.1. The results also demonstrate directly that little of the bolometric 
luminosity of the galaxy population arises from the faint end of the luminosity function, indicating a 
relatively flat faint-end slope of the IR luminosity function with a power law index of 1.2 ± 0.3. 

Subject headings: galaxies: evolution — galaxies: general — infrared: galaxies 


1. INTRODUCTION 

The rapid decay of the cosmic mean star formation 
rate (SFR) density from z = 1 to the present epoch has 
been convinci ngly e stablished over the last decade (see, 
e.g., TTonkius 120041 and references therein). The focus 
has now turned to characterizing the types of galaxies 
responsible for this decay. Rest-frame ultraviolet (UV) 
and optical emission-line studies indicate the decay since 
z ~ 1 is strongly influenced by the behavior of the rela- 
tively low-mass gala xies (e .g., Brinchma nn fc E llis et al. 
I200d Ju neau et al. Bundy et al. 120051 Bauer et 

al. I2OO5I: Wolf et al. T200^ . On the other hand, surveys 
in the thermal infrared (IR) have found many galaxies 
of intermediate and high mass with intens e, dee ply ob- 
scure d star formation (Fran ceschini et al. EoM Zheng 
et al. 12004 Bell et al. l2?)?)^ Perez-Gonzalez et al. 120051: 
LeFloc’h et al. 1200511 . However, the role of obscured 
star formation in lower-mass gjalaxie s remains unknown 
(although see Heavens et al. 12004 for a powerful and 
complementary approach). 

Among the ma ny observational SFR estimators (see 
Kennicutt Il998a| for a review), at high redshifts UV ra¬ 
diation is the most easily measured proxy for the SFR. 
It provides the basis for much of our understanding of 
the ev olutio n of the cosmic SFR density (e.g., Madau 
et al. mu Steidel et al. mni Sullivan et al. 
Hopkins et al. l20Mt Schiminovich et al. I200I . Yet, 
young stars are usually born in dust-rich environments; 
the dust absorbs the vast majority of the UV light and 
re-radiates this energy in the thermal IR. Therefore, to 
obtain a complete census of the bolometric luminosity 
from young stars, observations in both t he UV and IR 
are required (see, e.g., Gordon et al. mnai- However, the 
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available infrared facilities have lacked the sensitivity and 
imaging resolution to explore obscured star formation in 
low-luminosity, relatively low-mass galaxies, except for 
nearby examples. Thus, indirect estimates have had to 
be used for the IR outputs of such galaxies even at mod¬ 
erate redshift, using tools such as the empirical calibra¬ 
tion of UV color to extinc tion (e.g., Meurer et al. mn3 
Calzetti & Heckman ll99flll . or the locally-observed trend 
between Lir/Luv and luminosity (showing that more 
luminous galaxies tend to ha ve hig her extinctio n: e.g .. 
Wang & Heckman IT^ BelliOOl Buat et al. I2nn4i . 
These approaches have been applied to observations of 
high-redshift galaxies to derive the SFR per un it comov- 
ing v olume (e.g., Adelberger & Steidel l200rt Hopkins 

l2?ir)l . 

The goal of this paper is to measure directly the aver¬ 
age IR outputs of low-luminosity galaxies at intermedi¬ 
ate redshift. We use a deep 24/xm map of the Chandra 
Deep Field South (CDFS) from the Multib and Im aging 
Photometer on Spitzer (MIPS: Rieke et al. I2OOI . The 
24 fim map is limit ed by both photon and confusion noise 
(Dole et al. l2?inl . We find that stacking on the posi¬ 
tion of known intermediate redshift galaxies substantially 
reduces both sources of noise, allowing secure detection 
of average flux densities substantially below the conven¬ 
tional confusion limit. We have stacked 24/xm images of 
several thousand intermediate- and low-luminosity galax¬ 
ies at 0.1 < z < 1. We combine the resulting detections 
with the COMBO-17 data set to investigate the effects 
of dust obscuration in galaxies with mu < 24 mag out to 
z ~ 1. We find no evidence for evolution in the relation¬ 
ship between obscuration and the SFR over this range 
of redshift. Our results suggest that the IR luminosity 
function has a relatively flat slope (power law index of 
1.2 ± 0.3) toward low luminosities, at least out to z ~ 1. 

The paper is organized as follows. Section 2 describes 
the data and the selection of the galaxy sample. Section 
3 addresses the methods for stacking the images and de¬ 
scribes tests of their validity. Section 4 shows the results. 
Section 5 uses these results to refine estimates of the cos- 
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Fig. 1. — Illustration of Spitzer/MIPS 24 /rm images before and after PSF subtraction of individually-detected sources. Left an image 
section (10' X 10') in the CDFS; Right the corresponding residual image after removing all detected sources brighter than 83 /rjy (at the 
80% completeness). Two images are shown in the same greylevel scale. Sources seen in the residual images are fainter than 83/rJy. 


mic SFR density from z ^ 1 to the present epoch. Our 
work is summarized in Section 6. Throughout the pa¬ 
per we adopt a cosmology with i/p = 70kms“^ Mpc“^, 
Om = 0.3 and Oa = 0.7. 

2. THE DATA AND GALAXY SAMPLE 

In this paper, we combine ground-based optical data 
from the COMBO-17 photometric redshift survey with 
24 fj,m data from MIPS on board the Spitzer Space Tele¬ 
scope. As part of COMBO-17, the extended CDFS (cen¬ 
tered on 02000 = 03^32"'25*, i^aooo = -27°48'50") was im¬ 
aged in 5 broad and 12 me dium-b and optical filters from 
350 to 930 nm (Wolf et al. I2003D . T he CD FS catalog is 
available for public use (Wolf et al. I 200 I . The catalog 
presents astrometry and photometry for 63501 objects 
in an area of 30'.5 x 30'. Of these, photometric red- 
shifts and rest-frame luminosities in the standard John¬ 
son U,B,V bands and a synthetic UV band centered at 
280 nm"* are measured for ~9 000 classified galaxies with 
aperture magnitudes mu < 24 mag and z < 1.1. Red- 
shift uncertainties have been tested against large spectro¬ 
scopic datasets and are 5z/{l -I- z) ~ 0.02 at the median 
galaxy magnitude of R 2 2; astrometric uncertainties 
are ~ 0.1" (Wolf et al. [201)1 . 

MIPS observations of the CDFS at 24/rm were per¬ 
formed under MIPS Guaranteed Time Observations 
(CTOs). A rectangular field of 90' x 30' was observed 
in slow scan-map mode, with a total integration time of 
~ 1380 seconds per pixel. The 24 ^m image reduction 
wa s done with the MIPS Data Analysis Tool (Gordon et 
al. [200l . A final mosaic image was produced with reso¬ 
lution 1".25 per pixel and a Point-Spread Function (PSF) 
with Full Width at Half Maximum (FWHM) ~6". We 
used Sextractor (Bertin & Arnouts et al. 119961) for source 
detection. Because the typical size of an intermediate 
redshift galaxy is <C 6", we chose to fit all galaxies as un- 

The synthetic 280 nm passband has a square bandpass and 
40 nm FWHM. Rest-frame luminosities in the 280 nm bandpass 
are extrapolated using the best-fit galaxy templates for galaxies 
with z < 0.3. Rest-frame luminosities in U- and R-band are al¬ 
ways interpolated in the redshift range of interest, whereas V-band 
luminosities are extrapolations for galaxies with z ^0.7. 


resolved point sources. The few clearly extended sources 
were individually analyzed. An empirical PSF was built 
by stacking 18 bright point sources; this PSF was sub¬ 
sequently used for PSF fitting within 10" of all other 
objects^ The ALLSTAR routine in IRAF’s DAOPHOT 
was used to simultaneously fit multiple sources to the en¬ 
tire 24 fim image mosaic. Local background values were 
estimated to be the mode value of an annulus with an in¬ 
ner radius of 17" (5 FWHM) and width 20" centered on 
a given target®. The brightest sources were fit first and 
then removed. The fainter sources hidden by the bright¬ 
est sources were then detected and progressively included 
into ALLSTAR fitting. The flux of each detected object 
is derived from the PSF fits. We cut the 24 /im catalog 
at the 80% compl eteness level of 83 ^Jy (“5cr”, see Pa- 
povich et al. 120041 for the description of the completeness 
limit; the Papovich et al. catalog is very similar to that 
described here to within the errors). 

The MIPS and GOMBO-17 areas on the sky do not 
completely overlap. For the present analysis, we choose 
the subsample of galaxies contained within both the 
MIPS and the GOMBO-17 areas, avoiding the edges 
in the MIPS data. This sample covers an area of 
~800arcmin^ and contains 9785 objects with mu < 24 
(of which 7892 are galaxies in the redshift range 0.1 < 
z < 1). Before the core analysis of this paper, we con¬ 
structed a catalog of 24 /im detected sources that could 
be matched to GOMBO-17 sources with redshifts. We 
cross-correlated the two catalogs with a tolerance of 
1".25, or 3(7 positional uncertainty, as estimated from 
bright stars and compact sources. Of the 9785 GOMBO- 
17 objects at mu < 24, 1725 objects are MIPS 24/rm 
resolved sources. Of these, 1352 objects are selected in 
our sample as galaxies in the redshift range 0.1 < z <1. 

® Aperture corrections from this empirical PSF are similar to the 
updated PSFs in the “Spitzer Space Telescope Multiband Imaging 
Photometer for Spitzer (MIPS) Data Handbook” (version 2.3) to 
within the errors. 

® The use of median or mean values for the sky would make 
little difference for uncrowded sources; the mode is adopted as it 
is substantially more outlier-resistant than either mean or median 
in the case of sky annuli contaminated by sources. 














Stacking at 24 


3 



Fig. 2. — Reduction of the background RMS with increasing 
number of stacked images. Here fTt,g refers to the RMS of the flux 
within an aperture of radius 5”. For comparison, the dotted line 
shows the inverse square root of the stack number, normalized to 
the measured crbg at Ngtack = 10. 


3. MIPS 24/rM IMAGE STACKING 

In this section we describe the specifics of deriving the 
mean 24 fluxes for galaxy sub-sets through stacking. 
Basically, there are three steps: first, the removal of all 
individually-detected sources in the image; second, the 
stacking of the residual image postage stamps through 
averaging or medianing; finally, the addition of the in¬ 
dividual 24 ^m fluxes for all MIPS sources that coincide 
with the optical positions of the galaxy sub-set at hand. 

3.1. Stacking: method and testing 

Deep Spitzer surveys at 24 /rm, including these data, 
are limited by a combi nation of photon noise and confu¬ 
sion noise (Dole et al. l2nnl . The confusion noise at a 
given position on the sky is related to the unresolved ex- 
tragalactic population; increasing exposure time will not 
reduce this source of noise and therefore it sets a limit 
to one’s ability to individually detect galaxies in Spitzer 
24 /rm observations. However, the confusion noise does 
vary across the sky. If one is interested in estimates of 
the average flux for classes of objects (e.g. galaxies at a 
given redshift) with known coordinates, one can improve 
the detection threshold by stacking images centered on 
the object coordinates. The degree to which the confu¬ 
sion noise is thereby reduced depends on the degree to 
which the (confusing) sources are randomly distributed 
with respect to the population of interest. In the case 
of deep 24^m observations, the clustering is relatively 
weak (i.e., most close pairs of galaxies are projections 
rather than physical associations) and the background 
is already largely resolved (i.e., more than 50% of the 
total extragalactic 24/xm flux is in indivi dually -detected 
sources; Dole et al. l2nnl Papovich et al. 1200411 . making 
these data almost ideal for a stacking analysis. Hence, 
the combination of many images can substantially reduce 


the background noise and allow detection of the average 
flux of a given population, well below the canonical de¬ 
tection limit. 

How well this works in practice and how best to stack, 
we explore with a number of tests. As is clear from 
Fig.m a substantial portion of 24/im images have sig¬ 
nificant contributions from bright sources. In order to 
avoid that the “background” in stacked images is con¬ 
taminated by small numbers of bright sources, we first 
identified and PSF-subtracted all sources that were indi¬ 
vidually detected above the 83 ^Jy level. The right-hand 
panel of Fig.^shows a residual Image constructed in this 
way. We tested the efficacy of stacking these cleaned im¬ 
ages in two ways. First, we measured the RMS of the 
background in a randomly placed b" aperture (the same 
aperture used for the galaxy photometry) as a function 
of number of images stacked. Postage stamp images with 
a size of 1'.68 x 1'.68 were randomly extracted from the 
residual image and stacked. We quantify the background 
noise by modified, outlier-resistant measure of the vari¬ 
ance, defined as the half width of the range centered at 
the median value and including 68% of the sample values. 
Fig. El shows the decrease in crbg with increasing stack 
number: the effective background noise decreased as the 
inverse square root of the stack number. This behavior 
is as expected for ideal, uncorrelated random noise. The 
ideal behavior of the noise is a consequence of the flat 
number count spectrum at low flux levels at 24 /rm, and 
it has the benefit that a standard error analysis is useful 
on the stacked images. Only below Ubg ^ 0.4, which cor¬ 
responds to a “5cr” limit of ~3 to 4^Jy, does the effective 
background noise become significantly larger than that 
predicted by the inverse square law extrapolation (see 
also Fig. 0. Second, we tested the accuracy of recovering 
artificial sources from image stacks, by randomly placing 
large numbers of individually faint, identical sources into 
the source subtracted image. Fig. 0 shows the results: 
stacking is able to recover progressively fainter signals 
with larger stack numbers. Roughly speaking, for data 
of this quality (~ 1380 seconds per pixel slow scan map 
in a low background field), a 35 ^Jy signal can be recov¬ 
ered in a ^10-image stack with a flux accuracy of ±50% 
at the 3(7 detection level, and a 8 /iJy signal can be recov¬ 
ered with a ~100-image stack. In our stacking, postage 
stamp images are aligned to a few tenths of one pixel. 
Alignment to the centroid pixel, i.e., to an accuracy of 
one pixel, causes an additional 5% error on the final re¬ 
covered flux. 

3.2. Stacking of intermediate redshift galaxies 

To determine the characteristic thermal IR (24 /rm) 
flux for sub-sets of galaxies to z ^ 1, a sample was drawn 
from the COMBO-17 redshift survey, where it overlaps 
with the MIPS 24 ^m image. We divided the total sam¬ 
ple of 7892 galaxies into 9 equally-spaced redshift slices 
from 0.1 < z <1 and within each redshift slice, into sev¬ 
eral rest-frame H-band luminosity bins, scaled at each 
redshift to the characteristic absolute magnitude (Mg) 
of the lumin osity f unction derived by the VIMOS survey 
(Ilbert et al. l2nnl . Each luminosity bin covers one mag¬ 
nitude; adoption of a smaller magnitude bin would not 
change our results. We first illustrate which galaxies are 
individually detected through Fig. 0 which shows the 
rest-frame color U — V ot the sample galaxies versus the 
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Fig. 3.— Recovery of simulated sources with different input flux levels as a function of stack number. The intervals between stack 
numbers are identical in logarithm space. Assuming Poisson noise, the expected “5(7” limits for corresponding stack numbers are given in 
the left hand panels. The detection limit 83/xJy corresponds to stack number Agtack = 1- As the simulated sources in one stack set are 
identical, median stack and mean stack give almost the same recovery. Here median stack values are adopted. 


absolute B band magnitude. Open circles denote galax¬ 
ies individually detected above 83 /rJy by MIPS; dots de¬ 
note galaxies without individual detections. The dot¬ 
ted lines on the right-hand side of each panel show the 
value oi Mg in each redshift bin, while the dotted line 
on the left-hand side denotes the faintest magnitude of 
galaxies stacked in this paper. In the two lowest redshift 
slices, 0.1 < z <0.2 and 0.2 < z <0.3, the objects at 
the faint end have an abnormal color distribution, and 
have an unusually high incidence of detection s at 24 fj,m. 
Inspection of VVDS redshifts (LeFevre et al. HHoi and 
GEMS morphologies (Rix et al. 1200411 of galaxies in the 
low-luminosity, low-z bin argues strongly for a significant 
contamination by intrinsically more luminous galaxies at 
higher redshift. For very faint objects {mu ~ 23.5 mag) 
with a featureless blue continuum, it is difficult to deter¬ 
mine their redshift by either photometry or spectroscopy. 
Such objects are likely to be low- 2 : blue dwarf galaxies or 
luminous star-forming galaxies at high- 2 ; where 

the 4000A break is redshifted out of COMBO-17 broad 
/-band coverage). It is supported by the phenomenon 
that the high incidence of detections at 24 /rm occurs 


only in the low-luminosity and low -2 bins. Contamina¬ 
tion by high -2 objects to the low -2 objects is difficult to 
quantify. Accordingly, we excluded these faint objects 
from our investigation. It is worth noting that owing to 
COMBO-17’s /?-band selection, the faintest bin of the 
sample becomes incomplete at all redshifts 2 > 0.6; fu¬ 
ture work with a near IR-selected galaxy population can 
remedy this incompleteness. 

For each bin in magnitude and redshift, 1'.68 x 1'.68 
24/rm postage stamps from the PSF-subtracted image 
centered on the positions of the individually-undetected 
objects (see Fig. P) were cut out and then stacked. The 
size of the stamps is chosen to have sufficient area to 
properly estimate background. Two stacks were con¬ 
structed for each sub-sample: a stack in which each pixel 
is the mean, or alternatively the median of all contribut¬ 
ing pixels. The mean stack has the advantage that it 
includes the flux contributed from each galaxy, but has 
the disadvantage that it is more susceptible to residual 
flux from nearby sources. The median stack will under¬ 
estimate the total flux somewhat, as the contributions 
from bright, nearly-detected sources will not be incor- 
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Fig. 4.— COMBO-17 galaxy sample in the CDFS: each panel shows the rest-frame color U — V as a function of B band absolute 
magnitude for objects in 9 different redshift slices, increasing in redshift from top to bottom and left to right. In each panel, the right-hand 
side dotted line shows M* and the dotted line on the left shows the limit of the faintest bin considered in this work. The Mq adopted from 
Ilbert et al. is -20.0, -20.7, -21.1, -21.0 and -21.5 for redshift bins 0.1 < z <0.2, 0.2 < z <0.4, 0.4 < z <0.6, 0.6 < z <0.8 

and 0.8 < 2 < 1, respectively. Sources individually detected at 24 pm are shown as open circles; undetected sources are denoted by dots. 
The population of very faint 24 pm individually-detected sources at 2 < 0.3 consists of much more luminous galaxies at high redshift which 
were mistakenly assigned COMBO-17 photometric redshifts 2 < 0.3 (see text for more discussion). 


porated properly, but has the advantage that it is more 
robust to contamination from nearby sources. Further¬ 
more, the background values for the mean stack will be 
contaminated by sources, whereas the median will be 
more robust. Therefore, in the following we adopt two 
measures of the stack flux: the integrated 5" aperture 
flux of the median stack (using the median stack back¬ 
ground value for sky subtraction), a lower limit to the 
flux; and, the integrated central flux within an aperture 
of radius 5" of the mean-stacked image, using the back¬ 
ground estimate from the median-stacked image as the 
background, giving an upper limit to the flux. Finally, 


the empirically-derived correction of a factor of 1.881 was 
used to aperture correct the estimates of stack flux to to¬ 
tal. Fig. 13 shows 4 examples of median-stacked images. 
All median-stacked and mean-stacked images are elec¬ 
tronically provided in Appendix. Stack fluxes and their 
uncertainties are listed in Tabled along with total ob¬ 
ject numbers and stack numbers in each bin. Note that 
these mean flux estimates do not contain the flux contri¬ 
bution from the individual detections from each galaxy 
sub-sample. 

4. RESULTS 
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Fig. 5. — Examples of median-stacked 24pm images. The magnitude bin, redshift slice, integrated flux of the central target within an 
aperture of radius 5", /stack s-nd background fluctuation cri^g on the same aperture are labeled. Both /stack s-nd (Tt,g are given in units of 
pjy. No aperture correction is implemented. The size of each image is 1'.68 X 1'.68. In bottom-left image, the second Airy ring of the 
24 pm PSF is visible. 


4.1. Average SEDs and the contribution of 
individually-undetected sources 

Average luminosities at rest-frame 2800A, Johnson 
U, B, and 1^-bands and at observed-frame 24/tm were 
calculated for each (redshift, luminosity) bin of galax¬ 
ies. The 24 /im luminosity includes contributions from 
both individually-detected and individually-undetected 
sources. For the latter, the mean 24 /tm intensity was 
estimated by stacking. Bootstrapping was used to esti¬ 
mate the uncertainty in total luminosity, i.e. randomly 
extracting the same number of sources from the parent 
set, either individually-detected sources or individually- 
undetected sources, then summing the luminosity from 
each subset. Table |21 lists average monochromatic lumi¬ 
nosities (L^) in the rest-frame 2800 A, standard Johnson 
U,B,V, and observed MIPS 24 ^m bands for all mag¬ 
nitude bins in our investigation. Total UV luminosities 
derived from the 171/2800 A luminosities are given in unit 
of solar luminosity. Total IR luminosities and SFRs are 
also tabulated (see Sect. 14.21 & 14.81 for details). Fig. El 
shows the average 24 /rm luminosities in different R-band 
magnitude bins over the redshift range 0.1 < z < 1; lu- 
mino sities in some of the optical bands (from Wolf et al. 
IIdI are shown to delineate the overall SED; the ver¬ 
tical errorbar is 1 “ct” , which shows the 68% confidence 
region from bootstrapping. The horizonal bars show the 
band widths at each wavelength. Brighter galaxies in 
rest-frame B-band show higher mean 24 fj,m luminosities 


in all bands and redshifts than fainter galaxies. To place 
our measurements in context, three Spectral Energ y Dis- 
tribution (SED) templates from Devriendt et al. (Il999fl 
are arbitrarily scaled and shown alongside the observed 
SEDs. The three templates represent galaxies in three 
relevant modes of star formation: normal spiral galax¬ 
ies, starbursts (M82) and massive starbursts (Arp220). 
The importance of stacking in determining the average 
IR properties of distant galaxy populations can easily be 
discerned: the mean 24/im fluxes of the faintest galax¬ 
ies are an order of magnitude below the canonical “5cr” 
detection limit for individual detections for galaxies with 
z < 0.7, and are a factor of 5 below the individual de¬ 
tection limit for faint systems with 0.7 < z <1. Since 
the average 24 /im luminosity accounts for contributions 
from individually-detected sources, the mean 24 /im in¬ 
tensity of individually-undetected sources derived from 
stacking is indeed even lower. 

Fig.ia shows the 24 /im flux fractions derived from in¬ 
dividual detections and from stacking. For the bright¬ 
est galaxies (M* — 1 < Mb < M*), almost all of the 
flux arises from individual 24 /im detections across the 
redshift range 0.1 < z <1. The galaxies that are not 
individually detected have low average 24 /im luminosi¬ 
ties; this is expected given that a substantial fraction of 
the most luminous galaxies in rest-frame B-band are “red 
sequence” galaxi es wh ich contain little recent star forma¬ 
tion (Bell et al. l20f)l . For sub-B* galaxies, the major¬ 
ity of galaxies in each bin are not detected individually. 
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Fig. 6. — Average luminosities at rest-frame 2800A, Johnson U, B, F-bands, and at observed-frame 24^m for different magnitude bins 
(relative to an evolving Mg) and redshift slices. The label “—1 ~0” refers to the magnitude bin M* — 1 < Mb < M* -|- 0 and so forth. 
Three arbitrarily shifted SED templates are presented for comparison, taken from Devriendt et al. 419991). Horizontal dashed lines show 
the 83 /rjy detection limit. 


However, a large fraction of the total luminosity in each 
redshift and luminosity bin comes from the individually- 
undetected sources, especially for samples with z >0.4. 
For the faintest bin in each redshift slice, the contribution 
from individually-detected sources is typically small. 

COMBO-17 survey prese nts optically identified type 
I AGNs (Wolf et al. 120041) . Such objects are not in¬ 
cluded in our sample. However, optically obscured type 
H AGNs are not eliminated from our sample. Type H 
AGNs could potentially contribute to 24/im flux. To 
estimate the contamination from the type H AGNs, we 
cross-correlated our sample with the extended GDFS X- 
ray point source catalogs from IMs CPFS and 250 ks 
Chandra observations (Lehmer et al. l20?Tl . We iden¬ 
tified 123 galaxies in our sample as X-ray detected type 
H AGNs. We found that the average 24/rm luminosities 
decrease approximately 5% for galaxies at Mb < M*+2 
by removing the 123 objects. Low luminosity galaxies 
are almost not affected. For X-ray undetected type H 


AGNs (to 250 ks exposure) with optical counterparts at 
mu <24, radio and mid-IR investigations (Donley et al. 
1200^ Martinez-Sansigre et al. Hol) suggest that they 
are not more than X-ray detected ones. Then type H 
AGN’s contribution to the average 24 ^m luminosity is 
most likely up to ~10% in our analysis. This effect is 
negligible compared with uncertainties in estimating av¬ 
erage IR luminosities (typically ~0.3 dex; see Table |2Il. 
On the other hand, it is still unclear whether IR emission 
in obscured AGNs is mostly related to AGN activity or 
star formation in their host galaxies. We conclude that 
obscured AGNs have insignificant effect on our conclu¬ 
sions. 

4.2. The relationship between IR/UV and rest-frame 
B-band luminosity 

Local star-forming galaxies exhibit a clear correlation 
between the Lir/Luv ratio and rest-f rame R-ban d lumi¬ 
nosity (e.g. Wang & Heckman 1199(1 Bell 1200, in the 
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Fig. 7.— The relative contribution from individually-detected 
sources (heavily shaded region) at 24 fim compared to the contri¬ 
bution from undetected sources (lightly shaded region) for differ¬ 
ent magnitude bins in different redshift slices. In each panel, the 
left-hand side of shaded region shows the results for mean stack¬ 
ing scaled to unit and the right-hand side shows those for median 
stacking linearly normalized to mean-stack values, incorporating 
bootstrap error (the left-hand side shows the total mean-stack 
flux-hlcr bootstrapping error, and the right-hand side shows the 
total median-stack flux—Icr); the heavily shaded region shows the 
contribution from individually-detected sources (from the left-hand 
side to the right-hand side, the height of heavily shaded region rel¬ 
ative to overall shaded region shows minimal estimate of the con¬ 
tribution to maximal estimate). The total number of galaxies and 
number of individually-undetected galaxies stacked are shown in 
the upper and lower-right respectively. Dotted lines denote those 
magnitude bins having bias against red galaxies due to the i?-band 
selection effect. 


sense that lower-luminosity galaxies typically possess a 
lower Lir/Luv ratio albeit with large scatter. While we 
cannot test the validity of this relation galaxy-by-galaxy 
at intermediate redshift using the current datasets, we 
can investigate its validity of this relation in an average 
sense using stacking. 

Fig.ia shows the FirZ-Luv ratio as a function of rest- 
frame B-band luminosity at 0.1 < z <1. COMBO-17’s 
observed-frame i?-band selection allows access to galax¬ 
ies as faint as 0.04 L* (i.e., M* + 3.5) up to z ^0.45, 
0.1 L* (i.e., M* + 2.5) up to z ~0.65, and 0.25L* (i.e., 
M* + 1.5) up to z ~0.95. Owing to Spitzer’s lower sensi¬ 
tivity and larger PSF at wavelengths beyond 24 /rm, the 
vast majority of galaxies are undetected at longer wave¬ 
lengths. Therefore the total IR flux must be estimated 
from the 24 /xm flux by a rather large extrapolation, us¬ 
ing suites of local templates’^. Several sets of luminosity- 

^ We will tackle the technically-challenging stacking of longer- 
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Mg (M* + ) 

Fig. 8.— Relationship between the mean Fir/Fuv ratio and 
rest-frame B-band luminosity, scaled to Mg, covering 0.1 < z <1. 
The Bir/Luv ratio increases with luminosity, i.e., brighter objects 
show a higher Bir/Btjv ratio. 


dependent templates are used to estimate the total IR 
luminosity over Arest = 8 — 1000 /xm and its uncertainty 
(See Le Floc’h et al. 120051 for more details and an exten¬ 
sive discussion of observational tests o f 24 xx m-to-total-IR 
corr ection s: see also Cha ry fc Elbaz I 2 OO 1 I Papovich & 
Bell 1200^ and Dale et al. 120051 for more discussions). In 
Fig.lHI the shaded region shows the lower and upper lim¬ 
its for the Tir/Luv ratio mainly due to the uncertainty 
of Lir. Here the lower limit refers to median-stack total 
flux — 1(7 bootstrapping error — Icr IR luminosity error 
and the upper limit to mean-stack total flux -|- Icr boot¬ 
strapping error-h Icr IR luminosity error. The Tir/Tuv 
ratio is correlated with rest-frame R-band luminosity at 
all redshifts up to at least z ^ 0.8. Brighter galaxies on 
average show a higher Lir/Luv ratio, consistent with the 
sense of the local relation. As COMBO-17’s observed- 
frame i?-band selection may exclude red galaxies in the 

wavelength Spitzer data in an upcoming work. 
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-1 0 1 

log (SFR/Mp, yr^^) 

Fig. 9. — The mean Lir/Luv ratio as a function of SFR for galaxies in magnitude bins from M* — 1 < Mb ^ M* to M* + 3 < Mb < 
M* + 4. The heavy lines show increasing redshift, from left to right, with the minimum and maximum values of 2 at the corr espond ing 
ends of the lines. The local relation, recently re-derived using GALEX and IRAS data, is shown with dotted lines (Martin et ah I20Q5I) . 


faintest bin at z > 0.6, it is conceivable in these cases that 
the Lift/Luv ratio is underestimated (and especially at 
z > 0.8) due to a deficiency of red galaxies. 

4.3. The relationship between Lir/Luv and SFR at 

z < 1 

Hop kins et al. and Adelberger & Steidel 

Hoo) explore the relationship between UV-extinction 
and SFR, arguing for a universal correlation between ex¬ 
tinction and SFR. In this section, we explore this issue 
using our stacked dataset. 

For the data at hand, we estimate the SFR from a 
combination of UV (the directly-observed emission from 
young stars) and IR luminosities (dominated by repro¬ 
cessed UV light from young stars; e.g., Gordon et al. 
[ 2003 ). The IR luminosity in turn is esti mated from the 
24p,m luminosity. Following Bell et al. (1200511 . we esti¬ 


mate the “total” UV luminosity, from 1216 A to 3000 A, 
from COMBO-17’s 2800A monochromatic luminosity us¬ 
ing Lijy = 1.5 29 , 00 - This conversion assumes a stel¬ 

lar population with a constant SFR for 100 Myr and a 
Kroupa IMF. The same stellar population is then used 
to calibrate the total SFR: 

SFR/(M 0 yr-F = 9.8 x 10““ x (Lir-H 2.2Luv), (1) 

where Lir is the total I R lum inosity in units of solar 
luminosity (see Bell et al. uni for further description of 
the assumptions behind this SFR calibration). 

Fig. mi shows the Lir/Luv ratio as a function of SFR 
split into 5 bins of optical luminosity. The l ocal re lation- 
ship, recently re-derived by Martin et al. using 

GALEX and IRAS data, is also shown®. From the nearby 
universe back to z ~ 1, both Lir/Luv a-nd SFR increase 

® Martin et al. derive UV luminosities using vljp at 1500A, which 
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at a given optical luminosity relative to L*^. At all red- 
shifts and luminosities, galaxies follow approximately the 
local trend for Lir/Luv and SFR to within the (consider¬ 
able) uncertainties. In particular, as one considers fainter 
and fainter optical luminosity (from comparison of the 
different panels), galaxies have lower SFRs and Lir,/Luv 
values, again in agreement with locally-observed galaxies. 
Thus, to within the systematic uncertainties, we see no 
evidence against the proposition that 0.1 < z < 0.8 star¬ 
forming galaxies obey on average the local Lir/Luv“ 
SFR correlation. Beyond z ~ 0.8, the incompleteness of 
the sample makes such inferences inconclusive. This rela¬ 
tionship has been studied on a galaxy-by-galaxy basis at 
z ^ 0.7 by Bell et al. (2005) for individually-detected IR- 
luminous galaxies. These distant star-forming galaxies 
were found to approximately follow the local Lir/Luv“ 
SFR correlation, again to within the considerable un¬ 
certainties and with very large (> 0.5 dex) scatter. We 
stress that such scatter invalidates attempts at imple¬ 
menting the relationship to individual very luminous in¬ 
frared galaxies with SFR >2OM0yr“^. 

5. COSMIC IR LUMINOSITY DENSITY 

The cosmic IR luminosity density is essential to esti- 
ma te of t he cosmic star formation rate (e.g., Le Floc’h et 
al. HOQ^). It can be derived from the infrared luminosity 
function. Due to limited detection depth in the infrared, 
the faint end slope of the infrared luminosity function is 
poorly constrained, resulting in large uncertainty in esti- 
matin g the cosmic IR luminosity density (Perez-Gonzalez 
et al. 1200511 . Our stacked dataset probes flux levels be¬ 
low the individually-resolved flux in deep 24 ^m images. 
Hence, our results can provide new constraints on the 
cosmic IR luminosity density, which can be then used 
to constrain the slope of the faint end of the infrared 
luminosity function. 

We derived average IR luminosity for galaxies of known 
rest-frame H-band luminosity. By Comparing with B- 
band luminosity function, we are able to estimate the to¬ 
tal infrared luminosity density. Similar to Fig.|Hl we show 
the relation between Lir/z/Lj^ b ratio and the B— band 
magnitude in Fig. ^1 We fitted the behavior in this 
figure and used the fits to extrapolate to fainter lumi¬ 
nosities. The fits are given by the equations 

log Lm/vL^,B = -0.14 x Mr -L (0.6 ± 0.3) (2) 

and 

log LiYi /= —0.03 X ATb + (0.5 i 0.3) (3) 

for 0.1< z <0.4 and 0.4< z <0.7, respectively. Here 
the B-band magnitude is scaled to Mg. For redshift 
bins 0.7< z <1, we adopted the same slope as found for 
0.4< z <0.7 but left the normalization free, to derive: 

log Lui/vL^,B = -0.03 X Mb -h (0.8 ± 0.3). (4) 

The fits indicate that galaxies at the faint end are simi¬ 
lar to bright ones in the ratio, consistent with 

the finding of no correlation between the Btr / vL ,, r ra - 
tio and B-band luminosity by Wang & Heckman il996ll . 

yields values comparable to our definition of total UV for constant 
SFR populations. That work uses Lir = nm; these values 

are roughly a factor of two lower than our definition of IR lumi¬ 
nosity (see, e.g., Rell 120(1^ . We account for these differences in 
definition when comparing with the Martin et al. relation. 
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Fig. 10.— The relationship between mean Tir/Lb ratio and rest- 
frame B-band luminosity, scaled to M^, covering 0.1 < 2 : <1. The 
Tir/Tb ratio is roughly independent of luminosity to first order. 


Together with the re st-fra me ^-band luminosity func- 
tion of Ilbert et al. ll2005ll . we calculated the total IR 
luminosity density. In Fig. we show the total IR lu¬ 
minosity density as a function of redshift. Due to the 
observed-frame B-band selection, an increasing fraction 
of the total IR luminosity comes from the extrapolated 
part of the LF with increasing redshift. Since the B- 
band luminosity function is well dete rmine d in several 
deep redshift surveys (e.g. Faber et al. 1200511 . our results 
should not underestimate the cosmic infrared luminosity 
density. Similar results were obtained with the B-band 
lumin osity function of the DEEP2 survey (Willmer et al. 

Fig. E] exhibits that the cosmic luminosity den¬ 
sity increases by a factor of 9±3 from z = 0.1 to z = I. 

We compare our result with the total IR luminosity 
density given by the infrared luminosity function of Le 
Floc’h et al. II2005I1 in Fig. The two estimates of 
the total IR luminosity density are identical to within 
the errors oyer the redshift range 0.1< z <1. This is 
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0 0.2 0.4 0.6 0.8 1 

Redshift 


Fig. 11.— IR luminosity density as a function of red shif t. The thin solid line is the IR luminosity density attributed to galaxies whose 
population-averaged Tir/Tb ratio is known from Fig. [liJ By extrapolating the relation between Lir/Fb and Lb to the faint end, the 
total IR luminosity density derived from the B-band luminosity function is obtained and shown as the thick solid line. The thin dashed 
line shows the contribution from galaxies individually detec ted at 24 fizn (> 83 /^Jy). The thick dashed line shows the IR luminosity densit y 
given by the infrared luminosity function (Le Floc’h et al. |2003). The dot shows the local IR luminosity density (Takeuchi et al. |2003). 
The thick dashed line and thick solid line are well matched. 


important, as both estimates were derived in very dif¬ 
ferent ways. Le Floc’h et al. fit the evolving infrared 
luminosity function of individually-detected galaxies us¬ 
ing the same functional form as determined in the lo¬ 
cal Universe, incorporating luminosity and density evo¬ 
lution. In contrast, we directly detect the average flux 
of faint galaxies not individually detected by MIPS, and 
then use the trend in Lir/Lb as a function of Lb to esti¬ 
mate the IR luminosity residing in galaxies faintwards of 
our optical completeness limit (it should be noted that 
the shape of the optical luminosity function to deeper 
limits has been tested directly by the VVDS and e.g., 
in the HDF N and S, thus our extrapolation could be 
considered to be relatively well-posed). This comparison 
would tend to support Le Floc’h et al.’s claim that the 
faint end slope of IR luminosity function is unlikely to be 
very steep. We have tested this directly by re-fitt ing the 
IR luminosity functions of Le Floc’h et al. ll2005ll . using 


our estimated total IR luminosity (with our derived error 
bars) as a constraint. Such an analysis gives a power law 
slope of 1.2±0.3 (to be co mpared to the local value of 
1.23; Takeuchi et al. I2003t see also Perez-Gonzalez et al. 
120051 for a power law index of 1 to 1.3 up to z = 1 and 
beyond). 

6. DISCUSSION AND CONCLUSION 

Estimating the importance of dust-obscured star for¬ 
mation for intermediate-redshift normal, or even low- 
luminosity galaxies on a galaxy-by-galaxy basis is impos¬ 
sible through direct thermal IR observations with current 
technology. Using 24/xm images from Spitzer, we show 
that one can determine mean thermal IR fluxes and hence 
10-fold fainter average obscured star formation rates by 
stacking 24/xm images centered on the optical positions 
of known intermediate-redshift galaxies. We use a sam¬ 
ple of galaxies from the COMBO-17 photometric redshift 
survey of the Extended Chandra Deep Field South, which 
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provides astrometry, photometric redshift and rest-frame 
2800A, U, B, and F-band absolute magnitudes for thou¬ 
sands of 0 < 1 galaxies. We stack MIPS 24 /rm images for 
subsamples of galaxies in redshift slices and rest-frame 
B-band luminosity bins, allowing detection of average 
24 ^m fluxes down to <10/rJy, an order of magnitude 
deeper than those accessible on a galaxy-by-galaxy ba¬ 
sis. 

The mean total IR luminosity of these galaxy subsets 
is estimated from the observed 24 /rm luminosity, taking 
into account the uncertainty of IR SED shapes. Anal¬ 
ogous mean UV luminosities are derived, and average 
SFRs are estimated from the UV and IR luminosities of 
each subsample. We use these data to examine the corre¬ 
lations among optical luminosity, dust obscuration, SFR 
and redshift. 

We find that the correlation between dust obscuration, 
i.e. the ratio of Lir/LuV: and rest-frame B-band lumi¬ 
nosity seen in local star-forming galaxies holds over all 
redshifts z < I, with brighter galaxies showing a higher 
Uir/Auv ratio. Our averaged 24/rm detections show di¬ 
rectly that even in low luminosity galaxies (to 0.05B*) 
the majority of the bolometric luminosity from young 
stars is re-radiated in the thermal IR. Nonetheless, the 
decrease of Bir/Buv ratio with decreasing Br implies 
that globally star formation in faint objects is lower than 
the estimate one would derive the level of dust obscura¬ 
tion typical of normal galaxies. 

We explore the correlation between average Bir,/Buv 
ratios and SFRs for the different optically-selected sub¬ 
samples. Different subsamples populate different parts 
of the Bir/Buv“SFR plane; however, our data indicate 
that this correlation does not evolve much between z = 1 
and the present day. 

In closing, we briefly consider some of the factors de¬ 
termining the degree of dust obscuration indicated by 
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Uir/Uuv and roughly parameterized as an optical depth, 
r. The optical depth of a galaxy r oc Fgas Z a, where 
Ugas is gas density, Z is metallicity and a is a geometric 
term to account for the gas and dust distribution rela¬ 
tive to massive stars. Star formation intensity is strongly 
corre cted wi th gas density, e.g., the Schmidt law (Ken¬ 
nicutt IU)98H1 . Since our results are drawn from a large 
number of galaxies, geometric effects may cancel out, 
to first order. Then the correlation of Bir/Buvi or r 
with SFR (i.e., more intense star-forming environments 
show larger dust obscuration) may suggest that gas den¬ 
sity drives Bir/Buv correlation to a much greater extent 
than Z, at least over the magnitude and redshift ranges 
we consider in this work. 

Independent from the infrared luminosity function, we 
estimate the cosmic infrared luminosity density from the 
rest-frame B-band luminosity function. An increase by 
a factor of 9±3 is found for the comoving infrared lu¬ 
minosity density from z = 0.1 to z = l, consistent with 
Le Floc’h et al. il2nn^ . Based on our estimate of the 
comic infrared luminosity density, it is suggeste d that 
the infrared luminosity function of Le Floc’h et al. il2nnl 
is well determined at intermediate redshifts, supporting 
their claim that faint end slope of the luminosity func¬ 
tion is relatively flat. Our result suggests a power law 
index of 1.2 ± 0.3. 
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TABLE 1 

Total object numbers, stack object numbers, 24/jm stack fluxes and background RMS for 37 subsamples defined by 

REDSHIFT AND LUMINOSITY. 


z 

< Mb > 


stack 

Median-stack 

f stack 

(Jy) 

<^bg^ 

(Jy) 

Mean-stack 

/stack 

(Jy) 

(Jy) 

0.15 

-16.5 

206 

204 

3.37 

0.76 

2.82 

0.87 


-17.5 

134 

113 

7.33 

0.99 

7.13 

1.08 


-18.5 

93 

50 

11.90 

1.59 

12.82 

1.69 


-19.5 

69 

21 

12.00 

2.02 

14.26 

2.72 


-20.5 

17 

7 

22.80 

3.65 

18.73 

3.82 

0.25 

-17.2 

144 

141 

4.92 

0.97 

5.30 

1.02 


-18.2 

98 

88 

8.64 

1.10 

9.41 

1.31 


-19.2 

80 

51 

15.86 

1.37 

16.83 

1.76 


-20.2 

55 

20 

8.56 

1.89 

8.91 

1.99 


-21.2 

19 

12 

10.91 

3.22 

11.70 

3.42 

0.35 

-17.2 

175 

174 

1.02 

0.82 

0.77 

0.96 


-18.2 

165 

160 

5.47 

0.90 

5.37 

0.98 


-19.2 

116 

92 

12.89 

1.05 

13.49 

1.05 


-20.2 

57 

26 

17.42 

2.02 

17.03 

2.23 


-21.2 

20 

7 

6.34 

3.92 

4.55 

3.95 

0.45 

-17.6 

155 

153 

2.65 

0.79 

2.48 

0.73 


-18.6 

205 

192 

5.00 

0.68 

5.34 

0.80 


-19.6 

190 

165 

9.90 

0.75 

9.55 

0.83 


-20.6 

86 

50 

15.56 

1.39 

15.32 

1.34 


-21.6 

19 

11 

16.90 

2.71 

16.93 

3.44 

0.55 

-18.6 

362 

357 

3.71 

0.48 

3.69 

0.61 


-19.6 

315 

280 

10.02 

0.54 

9.85 

0.59 


-20.6 

200 

126 

14.53 

0.95 

15.37 

1.04 


-21.6 

67 

21 

12.20 

1.96 

11.78 

1.99 

0.65 

-18.5 

289 

282 

3.20 

0.67 

3.61 

0.77 


-19.5 

573 

533 

6.48 

0.44 

6.54 

0.50 


-20.5 

371 

263 

14.55 

0.65 

14.46 

0.66 


-21.5 

137 

64 

11.98 

1.27 

13.24 

1.29 

0.75 

-19.5 

576 

531 

6.68 

0.56 

6.98 

0.54 


-20.5 

343 

243 

14.78 

0.67 

14.80 

0.77 


-21.5 

130 

60 

10.94 

1.39 

12.52 

1.60 

0.85 

-20.0 

404 

324 

10.97 

0.53 

11.83 

0.61 


-21.0 

196 

111 

12.06 

0.92 

12.40 

1.02 


-22.0 

29 

13 

5.00 

2.38 

4.61 

2.55 

0.95 

-20.0 

327 

301 

8.44 

0.62 

8.74 

0.76 


-21.0 

250 

159 

12.25 

0.86 

12.22 

1.00 


-22.0 

39 

20 

6.49 

2.67 

6.75 

2.75 


'Within an aperture of radius . 







TABLE 2 

AVERAGE LUMINOSITIES AND SFRS FOR GALAXY SUBPOPULATIONS OF GIVEN MAGNITUDE BINS OVER THE REDSHIFT RANGE 0.1< Z <1. 
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z 

< Mg > 
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logLu 

log Lb 

log Ly 
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log Luv 
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log Lir 
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(W/Hz) 
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^Minimal estimate of IR luminosity and SFR from median-stack total flux — 1 bootstrapping error. 
^Maximal estimate of IR luminosity and SFR from mean-stack total flux -|- 1 cr bootstrapping error. 
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Stacking at 24 /im 

Appendix: MIPS 24 /im stack images 


15 


Here we provide MIPS 24 median-stack and mean-stack images of 37 subsets of galaxies in our sample. The 
images are available in the electronic edition of the Journal. 
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Fig. 12. — Median-stacked (left in a pair) and mean-stacked (right in a pair) MIPS 24 fj,m images for 37 stack bins in our sample. For each 
image, magnitude and redshift ranges are labeled at top-left; integrated central flux within aperture of radius 5^^ /stack? background 
RMS measured using the same aperture, crbg, are labeled at top-right. 
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Fig. 12.— Continued. 
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Fig. 12.— Continued. 
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